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ABSTRACT 

We present wide-field BVI photometry for about 11,500 stars in the low- 
met allicity cluster NGC 5466. We have detected the red giant branch bump 
for the first time, although it is at least 0.2 mag fainter than expected relative 
to the turnoff. The number of red giants (relative to main sequence turnoff 
stars) is in excellent agreement with stellar models from the Yonsei-Yale and 
Teramo groups, and slightly high compared to Victoria-Regina models. This 
adds to evidence that an abnormally large ratio of red giant to main-sequence 
stars is not correlated with cluster metallicity. We discuss theoretical predictions 
from different research groups and find that the inclusion or exclusion of helium 
diffusion and strong limit Coulomb interactions may be partly responsible. 

We also examine indicators of dynamical history: the mass function expo- 
nent and the blue straggler frequency. NGC 5466 has a very shallow mass func- 
tion, consistent with large mass loss and recently-discovered tidal tails. The 
blue straggler sample is significantly more centrally concentrated than the HB or 
RGB stars. We see no evidence of an upturn in the blue straggler frequency at 
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large distances from the center. Dynamical friction timescales indicate that the 
stragglers should be more concentrated if the cluster's present density structure 
has existed for most of its history. NGC 5466 also has an unusually low central 
density compared to clusters of similar luminosity. In spite of this, the specific 
frequency of blue stragglers that puts it right on the frequency - cluster My 
relation observed for other clusters. 

Subject headings: equation of state — globular clusters: individual (NGC 5466) 
- stars: blue stragglers — stars: evolution — stars: luminosity function 



1. Introduction 



The luminosity function (LF) is an observational tool used for analyzing the post- 
main sequence evolutionary phases of low-mass (~ 0.5-0.8M Q ) metal-poor stars in Galactic 
globular clusters (GGC). Because of their age and richness, GGC typically contain hun- 
dreds of stars that have evolved off the main sequence. The numbers of stars in evolved 
phases are directly related to the evolutionary timescales and fuel consumed in each phase 
( Renzini &: Fusi Peccilll988l ). so that they present us with an opportunity to test this aspect 
of stellar evolution models. 

The results of the most st r ingent tests have been mixed. Repeated stu d ies of the metal- 



poor cluster M30 (IBoltd 11994 : iBergbuschl Il996l ; iGuhathakurta et al.l Il998l ; ISandquist et al. 



19991 ) have f ound an e xcess number of red giant branch (RGB) stars relative to main sequence 
(MS) stars. iStetsonl (119911 ) also uncovered an apparent excess of stars in a combined LF of 
the metal-poor clusters M68, NGC 6397, and M92. However, the LFs of more metal-rich 
clusters show no discrepancy (M5: Sandquist et al. 1996; M3: Rood et al. 1999; M12: 



Hargis, Sandquist, & Bolte 2004). In a survey of 18 clusters, IZoccali &; Piottol (120001 ) found 
good agreement with model predictions with the possible exception of clusters at the high 
metallicity end. 

In this paper we present BVI photometry of NGC 5466, a high galactic latitude globular 
cluster (/ = 42.2° and b = 73.6°), located in the constellation of Bootes (a = 14 h 05 m 27?4, 
5 = +28°32'04" at a distance of R=15.9 kpc; Harris 1996). NGC 5466 is a loose cluster 
(r c = l'.6A) with extremely low metalli city ([Fe/H]= —2.22) and subject to little or no 
reddening, (E(B - V) ~ 0) jHarrislll996h . 



In §2, we describe the process leading to the calibrated photometry, and compare with 
previous studies of the cluster. In §3, we compare the observed color- magnitude diagram 
and observed luminosity function with theoretical models, focusing on the relative number 
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of stars on the lower RGB and around the MS turnoff. Finally, in §H we present a new 
examination of the blue straggler population of NGC 5466. 



2. Observations and Data Reduction 

The data used in this study were obtained with the Kitt Peak National Observatory 
(KPNO) 0.9 m telescope (0"68 pix" 1 ) on the nights of UT dates 1995 May 4, May 5, and 
May 9. A complete list of the image frames, exposure times, and observing conditions is 
given in Table [U 

The images obtained on the three nights were processed using tasks and pack- 

ages. The reduction involved subtraction and trimming of the overscan region of all images, 
subtraction of a master bias frame from flats and object frames, and flat fielding of the 
object frames using images taken at twilight. Profi le-fitting; photometry was done using the 



DAOPHOTII/ALLSTAR programs (1Stetsonlll987f ). 



We also reduced archival ground-based photometry of the cluster core taken with the 
High- Resolution Camera (HRCam) on the 3.6 m Canada- France-Hawaii Telescope (CFHT). 
The V and / images were taken 30 and 31 May 1992 (observers J. Heasley and C. Christian), 
and have not previously been described in the literature. The CCD images had 1024 x 1024 
pixels, 0'.'13 per pixel, and excellent seeing (0.4 - 0.5 arcsec). The images were reduced using 
the archived bias and twilight flat frames, and following a procedure similar to that for the 
KPNO data. This allowed us to get excellent photometry to 2 magnitudes below the turnoff 
in the cluster core. These images were used entirely for blue straggler identification (see §1]). 



2.1. Calibration against Primary Standard Stars 



The conditions at KPNO on 1995 May 9 were photometric, and lLandoltl standard star 
fields were observed at a range of air masses to determine photometric transformation coeffi- 
ci ents. Th e stand ard values used for the calibration were chosen from the large co mpilation 



of IStetsonl (120001 ) . which is set to be on the same photometric scale as the earlier lLandolt 
Jl992|) values. 

We conducted photometry on the standard stars and isolated cluster stars using multiple 
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synthetic apertures. We then used the DAOGROW ( Stetson! 1990 ) program to construct 
growth curves to extrapolate measurements to a common aperture size. Using the CCDSTD 
program, the standard star transformation equations were found to be: 

b = B + a + (-0.069 ± 0.005)(£ - V) + (0.255 ± 0.014)(X - 1.0) 
v = V + b + (0.027 ± 0.004X5 -V) + (0.172 ± 0.010)(X - 1.0) 
i = I + c + (-0.013 ± 0.005)(V -T) + (0.145 ± 0.018)(X - 1.0) 

where X is the airmass, v, b and i are instrumental magnitudes, and V, B, and I are standard 
mag nitudes. Thes e calib ration equ ations are dif f erent than those used for our analysis of 
M10 iPollard et al.l (120051 ) and M12 lHargis et al.l (120041 ). which were observed on the same 
night, because our /—band exposures did not go as deep as the B and V exposures. As a 
result, the (B — V) color was a better choice for calibrating the V photometry down to our 
faintest observed stars. The calibrated measurements for the standard stars are compared 
with catalogue values in Fig. [TJ 

We note that there is slight evidence of trend in the residuals for the I band versus 
magnitude, which might indicate nonlinearity. This impression is caused by one observation 
of the PG1323-086 field. We did, however, have an additional observation of the same field 
on the same night having the same exposure time that does not show the same (small) trend. 
Because we do not have any reason to eliminate the frame and because its elimination has a 
minimal effect on the transformation coefficients, we have decided to retain the measurements 
from the image. 



2.2. Calibration against Secondary Standard Stars 

Aperture photometry for 165 cluster stars was used to calibrate the point-spread function 
(PSF) photometry for the cluster. These secondary standard stars were chosen based on 
relatively low measurement errors and location in relatively uncrowded regions of the cluster. 
They were chosen from the asymptotic giant branch (AGB), upper RGB, and horizontal 
branch (HB) of the cluster in order to cover the entire range of colors covered by cluster 
stars. 

The PSF-fitting photometry for the three nights of data was combined and averaged 
after zero-point differences among the frames had been determined and corrected. The zero- 
point corrections to the standard system were determined after fixing the color- dependent 
terms at the values measured in the primary standard star calibration. (This was also done 
in our studies of M10 and M12.) In Fig. [21 it can be seen that this procedure does not 
introduce systematic color- or magnitude-dependent errors. 
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2.3. Comparison with Previous Studies 



We co mpared our ph otometric data set to those of IJeon et al.l (120041 ). iRosenberg et al. 



(120001 ) . and IStetson (120001) . The mag nitude and c olor compar isons (BVI in this study versus 
VI in IStetsonl and IRosenberg et al.l . and BV in IJeon et all ) as function of magnitude and 
color are shown i n Figs . [3]- El Though our calibrated magnitudes are slightly brighter than 
those of IStetsonl (120001), the diffe r ences are small, and there is no color trend. The offsets 
compared to the IRosenberg et all (120001) are lar ger, but again there are no clear color trends. 



The offsets compared to the 
are slight trends with color. 



Jeon et all (120041 ) data are also significant, but more notable 



2.4. Calculation of the Luminosity Function 



Artificial star tests were performed to empirically measure the precision of our photom- 
etry and to correct for incomple teness in the detection of stars. We followed the procedure 
described by lHargis et all (120041 ) for the calculation of incompleteness corrections as function 
of position and magnitude. 

The inputs used for producing the artificial star tests were the reduced B and V 
CC D frames, PSFs for ea ch object frame, fiducial lines, and an estimate for the initial 
LF (jSandquist et allll996l ). Artificial stars were randomly placed in cells on a spatial grid 
and the entire grid was t hen shifted randomly fr om run to run in order to ensure the whole 
imaged field was tested (IPiotto k, Zoccalil 119991 ) . Each star was placed in a consistent po- 
sition relative to the cluster center on each image. If a detected star was found to coincide 
with the input position of an artificial star, it was added to the archive. The new images 
were reduced using the same procedure applied to the original data set. In this study, a 
total of 100,000 artificial stars from 50 separate runs were added. The number of artificial 
stars per trial was chosen so that the effects of crowding on the photometry was qualitatively 
unchanged. 

The recovered artificial stars were used to calculate 1) median magnitude and color 
biases (5y and 5b-v, where 5 = median [owtpwt — input]), 2) median external error esti- 
mates (a ext {V) = median[5y — median(<5y)]/0.6745 and o ext (B — V)), and 3) total recovery 
probabilities (F(V), which is the fraction of the stars that were recovered with any output 
magnitude) in bins according to projected radius and magnitude. The values for the above 
quantities are plotted in FigsJ2]-[EJ 



Finally, an initial estimate of the "true" LF and the error distribution, magnitude biases 
and the total recovery probability (F) were used to compute the completeness fraction / 
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(the ratio of the predicted number of stars to the actual number of observed stars). The 
completeness fraction results are shown in Figure [9j We then interpolated to compute / 
for the radial distance and magnitude of each detected star. For each observed star f^ 1 
was added to the appropriate magnitude bin to determined the observed LF. (Note that the 
completeness fraction was set to 1.0 for star brighter than the turnoff.) The observed LF 
along with the upper and lower 1 a error bars on log N are listed in Table [2j 



3. Discussion 



3.1. Reddening, Metallicity, Distance Modulus, and Age 

Bec ause NGC 5466 resi des at high galactic latitude, it suffers little if any reddening. 
Though ISchlegel et al.l (119981 ) found a redd ening of E(B — V) = .02 from the maps of dust 
IR emission, we adopted E(B — V) = 0.0 ([Rosenberg et al.ll2000l ). For our interests in this 
paper, the small difference is of small importance. Most of the comparisons below between 
observations and theory are relative, in which reference points (like the turnoff) are used to 
determine magnitude and/or color shifts. This has the benefit of minimizing the influence 
of uncertainties in reddening and distance modulus (see below). 

As for abundances, there is only one high-resolution measure ment for a cluster star, and 
it is for the anoma l ous C epheid V19. iMcCarthv fc Nemed Jl997h find [Fe/H]= -1.92 ±0.05, 
while iPritzl et al.l (120051 ) find [Fe/H] = —2-05 u sing the same data. Typically quoted metal- 
licity values include [Fe/H] = —2.17 (jZinn!ll980l ) from photoelectric phot ometry of integrated 
light in selected filter bands, and [Fe/H] = —2.22, which was derived b y lZinn fc Westl (119841 ) 
(ultimately from low- resolution s pectral scans bylSearle fc Zinnlll978l ). When converted to 
the widely-used metallicity scale of lCarretta fc Grattonl (119971 ). this becomes [Fe/H] = —2.14. 
More work could certainly be done on the composition of NGC 5466 stars, but the evidence 
so far points to an abundance [Fe/H] < —2.0. Though the range in the above quoted metal- 
licity values is relatively large for a globular cluster, the exact value is not critical for our 
purposes since we will primarily be concerned with relative comparisons. 

Our photometry does not extend f aint enough to derive a new distance modulus from 
subdwarf fitting to the main sequence. iHarrisl (119961 ) obtained (m — M)y = 16.0 by cali- 
brating the observed luminosity level of the horizontal branch with the relation My(HB) = 
0.15[Fe/H] + 80 an d adopting a reddening, E{B -V)=0.0 and a metallicity, [Fe/H] =-2.22. 
Ferraro et al.l (119991 ) determined distance moduli (m — M )y = 16.16 from their zero-a ge HB 
estimate, assuming no reddening and metallicity on the ICarretta fc Grattonl (119971 ) scale. 
We will consider distance moduli in this range. 
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Most previous age estimates of NGC 5466's age have it older than the recent deter- 
mination of the age of the univer se (13. 71°' 17 Gyi) obtained by the Wilkinson Microwave 
Anisotropy Probe (WMAP) team (jSpergel et al.ll2006l ). Recent homog eneous studies of GGC 



indicate that NGC 54 66 is coeval with clusters of similar metallicity (ISalaris Weiss! 12002 



Rosenberg et al.ll2000l ). As a result, we will primarily consider ages in the range of 12 to 13 
Gyr. 



3.2. The Color-Magnitude Diagram 



The color-magnitude diagrams (CMDs) for NGC 5466 show well-defined RGB, AGB, 
and HB sequences (see Fig. [10]), and stars extending from the tip of the RGB down to 
V ~ 22.5. Fiducial sequences for the MS and lower RGB were determined from the mode of 
the color distribution of stars in magnitude bins. The SGB position was determined using 
the magnitude distribution of the stars in color bins. The fiducial line for the rest of the 
RGB was obtained from the mean color of stars in magnitude bins. The fiducial points are 
listed in Table [31 

A comparison of the fiducial p oints derived for NG C 5466 with theore tical isochrones for 



a ran ge of ages from the Teramo ( Cassisi et a. 



20061 ) , and Yonsei-Yale (IDemarque et al.l 12004 ) groups is displayed in Figures [IT] and 



20041 ). Victoria- Regina (IVandenBerg et al. 



The isochrones have been shifted in color and magnitude (aligning the turnoff colors and 
the magnitudes of the main sequence point 0.05 mag redder than the turnoff) according 
to the technique of IVandenberg et al.l (jl990l ). This has the advantage of removing some of 
systematic uncertainties associated with the color-T e jj transformations. [In our comparisons, 
we found that the Yonsei-Yale models could not match the fiduc i al line for any reasonable set 
of input parameters when the transformations of lLejeune et al.l (119981 ) were used. Therefore, 
we only utilize models using the Green, Demarque, & King (1987) transformations. Even 
then, we could not find a match with the slope of the upper giant branch for reasonable 
metallicities ([Fe/H] < —1.9).] On the whole, the shape of the fiducial matches the models 
well on the main sequence, subgiant branch, and lower giant branch. Neither the Teramo 
nor the Victoria-Regina models include element diffusion processes, while the Yonsei-Yale 
isochrones only include He diffusion. However, differences in T e //-color transformations are 
likely to be the cause of some of the differences seen. 
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3.3. The Luminosity Function 

The number of stars at a given luminosity in post-main-sequence phases is directly pro- 
portional to the lifetime spent at that luminosity. It is well known that the LF of the RGB 
probes the chemical stratification inside a star because the hydrogen abundance being sam- 
pled by the t hin hydrogen-fusion shell a ffects the rate of evolution, and hence star counts 



on the RGB (IRenzini fc Fusi Peccil Il988l ) . Setting aside the short pause at the RGB bump, 
RGB evolution accelerates in a very regular way that is ultimately related to the structure 
of degenerate core. The relationship between core mass and radius forces the fusion shell to 
function at strictly controlled density and temperature conditions, which leads to a relation- 
ship between core mass and luminosity. This causes the LF to be particularly sensitive to 
certain physical details, which we discuss in §3.3.11 below. 

Fig. [13] shows the observed luminosity function compared to theoretical LFs for the la- 
beled values of metallicity, exponent of the initial mass function, and a range of age estimates 
for NGC 5466, assuming our preferred distance modulus (m — M)y = 16.00. The theoret- 
ical models were normalized to the observed LF at V ~ 21.3 (sufficiently faint that stellar 
evolution effects are minimized). The models agree well with the observed LF, implying an 
age of approximately 12 - 13 Gyr. 



3.3.1. Relative RGB and MS Numbers 



Gallart et al.l (120051 ) recently discussed theoretical luminosity functions calculated by 



different groups. One of the primary differences they noted was in the number of giant stars 
relative to main sequence star s. In order to show t hese d ifferences in a parameter- independent 
way, we follow the method of IVandenberg et al.l (Il990l ). In Fig. [T4], the theoretical LFs were 
shifted so that a point on the main sequence 0.05 mag redder than the turnoff color point 
matched the corresponding point on the cluster fiducial line. The reason for using the point 
(Vro+o.os) rather than the MSTO itself is that the MS has a significant slope and curvature 
at this point, m aking it possible to accu rately measure the point in both observational data 
and isochrones (IVandenberg et al.lll990l ). The theoretical models were normalized to the two 
bins in the observed LF on either side of the turnoff (V = 19.83 and 20.13). As can be seen, 
age-related differen ces between the theoretical LFs nearl y disappear when this procedure 
is applied (see also Stetson 1991 ; Vandenberg et al. 19981 ). However, when different sets of 
models are compared, there are small differences in the number of RGB stars relative to 
MS stars, with the Victoria-Regina models predicting the smallest number of giants and the 
Yonsei-Yale models predicting the largest. 
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To quantify the differences, we computed the ratio of t he number of stars o n the lower 
giant branch to the number of stars near the main sequence. iPollard et al.l (120051 ) introduced 
this ratio and showed that it is insensitive to age and heavy-element abundance. For the 
main sequence population, we used star counts in the two bins on either side of the turnoff 
(19.682 < V < 20.282), and for the red giant branch we used the counts in the range 
16.982 < V < 18.482. We derived model values from the same magnitude ranges relative to 
the Vto+o.05 point on the main sequence. Values are compared in Table HI The error in the 
observed value is dominated by Poisson statistical scatter. The Yonsei-Yale models are in 
best agreement with the observations, the Victoria-Regina models are out of agreement by 
more than 2 er, and the Teramo models are in between. 

It is worth examining the possible causes of this difference both because it may help 
improve the physics inputs for the models and because red gi ant stars are some of the largest 
contributors to the integrated light of old stellar populations. iGallart et al.l (120051 ) t abulated 
most of the main physics inpu ts for the most widely-used model set^]. Earlier studies ( Stetson 
199ll ; IVandenberg et al.lll998l ) have shown that the LF-shifting method used above eliminates 
nearly all sensitivity to model input parameters like mass function, age, convective mixing 
length, and composition inputs with the exception of helium abundance, which we examine 
first. 



Older models (e.g. Fig. 9 of Ratcliff 1987, Fig. 7 of Stetson 1991, Fig. 3 of lVandenberg et al 



19981 ) seem to agree that an increase in initial helium abundance Y in non-diffusive models 
by 0.1 results in an increase in the relative number of stars on the RGB (more precisely, a 
reduction in the relative number of main sequence stars) by about 0.07-0.08 in logiV. The 
Teramo models (Y = 0.245) predict about 12% more giant stars relative to main sequence 
stars compared to the Victoria-Reg ina models (Y = 0.235)|. This difference corresponds 
to a shift of 0.05 in logiV, which is about an order of magnitude too large for the helium 
abundance difference. 

The Yonsei-Yale models have the lowest assumed helium abundance (Y = 0.23), but 
are the only set of the three that include helium diffusion. The inclusion of helium dif- 
fusion reduces the age derived from the turnoff of a globular cluster by about 10-15% 



2 Since the lGallart et al.l ( 20051 ) review was published, the Teramo group found an error in the evolution 
scheme they used on the giant branch, which brings their models into better agreement with other groups. 
We use their updated models in the comparisons here 

3 The Teramo models also assume a larger a-element enhancement ([a/Fe] = 0.4) than the Victoria-Regina 
models ([a/Fe] = 0.3), which would tend to reduce the number of giants relative to main sequence stars. 
However, because the relative number of RGB and MS stars is not sensitive to small differences in heavy 
element abundance, this difference is probably unimportant. 
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( jProffitt fc Vandenbergl Il99ll ; IStraniero et al.l 119971 ; IVandenBerg et al.l 120021). thanks to the 
inwar d motion of helium. According to theoretical models (e.g. Fig. 8 of lProffitt fc Vandenberg 
199ll ). diffusion has a small effect on the LFs (~ a few times 10~ 2 in logiV, increasing for 
increasing age), but it does increase the number of giants relative to MS stars. He diffusion 
reduces the total core hydrogen fuel supply available to an MS star, but in itself this does not 
strongly modify the LF, just changes the brightness of the turnoff. This magnitude change 
is eliminated in our LF shifting procedure. The chemical composition profile left in the star 
after it leaves the main sequence has a greater impact. Diffusion reduces the H abundance in 
th e fusion regions, thereby d ecreasing the evolutionary timescale. According to the models 
of lProffitt fc Michaudl (119911 ). the changes to the He abundance profile are most considerable 
immediately below the surface convection zone, and just outside the nuclear fusion regions 
(where the composition gradient slows the inward settling of helium). However, over most 
of the star, the changes in Y are limited to 0.01 - 0.02. Because the core portion of the 
composition profile is consumed on the subgiant branch, the evolution timescales for giant 
stars are only affected in a minor way, and the appearance of a deep convection zone on 
the lower giant branch wipes out most of the effects of diffusion for the upper giant branch 
evolution. 

In spite of this, the Yonsei-Yale models have almost 23% more giants than the Victoria- 
Regina models, and more than 9% more giants than the Teramo models (relative to MS 
stars). The lower helium abundance in the Yonsei-Yale models compared to the other models 
should partially counteract what effects helium diffusion might have had on the RGB/MS 
ratio as well. Thus, it appears that neither He abundance nor He diffusion can completely 
explain the differences between the Yonsei-Yale models and the other groups. 



We can ask whether the LFs show similar disagreements at other metallicities. lHargis et al 



(120041 ) made comparisons between the Victoria-Regina a nd Yo nsei-Yale theoretical L Fs and 
observational LFs for the clusters M3 faood et al.lll999h . M5 JSandquist et alill996h . M12, 
and M30. The overall impression from those comparisons was again that the Yonsei-Yale 
models (having He diffusion) predict more giant stars relative to main-sequence stars than 
do the Victoria-Regina models. In Fig. [151 we compare the LFs for these clusters with 
the Teramo models. The degree of agreement or disagreement can be quantified with num- 
ber ratios of lower RGB and MSTO stars, similar to the ones we computed earlier for 
comparisons with NGC 5466. Our calculations are shown in Table HI As can be seen, uncer- 
tain ties in the metallicity scale have some effect on the comparisons w ith the observations . 
The ICarretta fc Grattonl (119971 ) scale has higher [Fe/H] values than the lZinn fc Westl (119841 ) 
scale, and thus results in lower number ratios. 



Fig. [H] shows the results of comparing the observed ratios with the models for different 
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[Fe/H] scales. On the IZinn fc Westlll984l scale (right panels), the observed values seem to 
be in agreement to within about 1 — 1.5cr for most of the models, with the exception of 
the the lowest metallicity clusters (M30 and NGC 5466) and the lowest helium abundance 
( Victoria- Regina) models. On the ICarretta fc Grattonl 119971 scale, the Yonsei-Yale models 
have the best overall agreement, although the Teramo models only deviate noticeably for 
the lowest metallicity clusters. The Victoria-Regina models predict too few giants for all of 
the clusters. 

The differences from model to model (as opposed to models versus observation) point 
toward deficiencies elsewhere in the physics or computational algorithms used in the stel- 
lar evolution codes. @ The RGB LF is a robust prediction of the models because there is 
a strong core mass — luminosity relationship: the conditions in the hydrogen fusion shell 
of the giant are strongly dependent on the structure of the degenerate core and are al- 
most independent of the details of the mass or structure of the envelope. As a result of 
this, we can focus on facto rs affecting core structure. (As a non-standard physics exam- 
ple, IVandenberg et al.lll998l describe the way in which core rotation relieves a giant star of 
some of the need to support itself by gas pressure, which reduces the core temperature and 
lengthens the evolutionary timescale.) Because model-to-model differences appear even on 
the faint end of the giant branch, we can set aside factors that only become important to the 
structure of t he star near the tip of th e giant branch [such as neutrino losses and conductive 
opacities; see lBjork fc Chaboyerll2006l . for example], even though there are significant theo- 
retical uncertainties in these quantities. Nuclear reaction rates in the fusion shell can also 
be neglected, par t ly bec ause the uncertainties in the reactions appear to be relatively small 
Adelberger et aJj (119981 ). but also because small changes in the reaction rates require only 
tiny changes in the shell temperature to get the same energy production. This leads us to 
examine the equation of state (EoS) in the core. 

Although the behavior of degenerate electrons is thought to be very well understood, 
their interactions with nuclei can have a measureable effect on the pressure. Particles of like 
charge tend to cluster together, which modifie s the free energy o f the g as and reduces the gas 
pressure for given density and temperature. lHarpaz fc Kovetzl (119881 ) looked at the effects 
of the inclu sion of Coulomb int eractions on giant stars, and their results are corroborated 
by those of ICassisi et al.l (120031 ) . They found that for a given core mass the fusion shell 
temperature was higher when the Coulomb interactions were included, which leads to faster 
processing of hydrogen. Thus this is another example (like core rotation) where modification 
of the pressure support of the core affects the evolutionary timescales, which results in 



4 Fo r a comparison of these physical inputs for the different theoretical groups, see Table 1 of 
d2005h . 



Gallart et al. 
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changes to the luminosity function. The Coulomb corrections to the pressure become more 
important with increasing density for the core, but are small compared to the contribution 
of the degenerate electrons. 

All of the model sets we have considered here incorporate Coulomb interactions in some 
form. The Teramo group used the most sophisticated "EOS1" version of the FreeEO£@, 
which incorporates Coulomb corrections in a form that matches limits in both the weak 
(Debye-Hiickel) and strong (once-component plasma) Coulomb interaction limits as well 
as (less importantly) electron exchange interactions. The strong interaction limit is most 
relevant for giant star cores since the strong interaction parameter 

C 2 e 2 

r = - > 1 

r kl 

where ( is the rms nuclear charge and ro is the average internuclear distance. 



The Yonsei-Yale group use d the OPAL EoS tab les ([Rogers et al.lll996l ). but falls back 
on the group's older EoS [e.g. (IGuenther et al.lll992l )] for conditions for high densities and 
temperatures outside the OPAL tables. While the most recent OPAL tables probably contain 
the most complete physical description of the Coulomb effect, the OPAL tables they used 
(Y.-C. Kim, private com munication) were computed prior to recent improvements to account 
for relativistic electrons (IRogers &: Nayfonovl 120021) . and as a result cut off at logp > 5.0. 
The Yale EoS at higher densities only includes the Coulomb effect in the weak Debye-Hiickel 
limit, which for the highest densitie s in the core. The Victoria- Regina models also use 
a modified version of the EFF EoS (lEggleton et al.lll973T). with a c orrection for Coulomb 
corrections in the weak Debye-Hiickel limit (jVandenBerg et al.ll2000l ). 



The differences in the implementation of the Coulomb effect may explain the fact that 
the Teramo models generally predict more giants (relative to the main sequence) than the 
Victoria-Regina models do. However, the smaller Coulomb corrections in the Yonsei-Yale 
models would tend to result in fewer giants than the Teramo models (although the effects of 
helium diffusion work in the opposite direction). So, we are unable to completely reconcile 
the differences in the luminosity functions from the three groups. 

Obviously more detailed study is needed by all of the modelling groups to identify the 
causes of the differences, but such a study is beyond the scope of this paper. Still, we 
believe that helium diffusion and strong interaction Coulomb corrections are physical effects 
that should be considered first. There is, for example, good evidence from helioseismology 



3 FreeEOS is available at http://freeeos.sourceforge.net/ and the discussion of the implementation of the 



Coulomb effect can be found at |http://freeeos.sourceforge.net/coulomb.pd4 
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for helium diffusion i n the Sun iBahcall et al.l (119951 ). despite the surface convection and 



meridional flow (e.g., lHathawayl Il996l ) . It is expected that helium diffusion should also act 



in globular cluster stars. A detailed study of the effect of eq uation of state uncertainties 
has yet to be done [see, for example, iBjork fc Chaboyerl (120061 ) for a study of uncertainties 
in other physical inputs]. Use of FreeEoS would make a study of equation of state effects 
most stringent since it appears to be capable of modelling the most sophisticated tabular 
EoS (OPAL), while also having the flexibility to allow individual bits of the physics to be 
"turned off". 

As a final warning about the observations, we should remember the LF of the clus- 
ter M10. iPollard et al.l (120051 ) found that unusual variations in numbers of RGB stars at 
different brightness levels in M10 (a virtual twin to M12). In particular there seemed to 
a significant excess in the number of stars near the RGB bump in brightness, while the 
lower RGB appeared normal (compared to Victoria-Re gina and Yonsei- Yale models). A 
similar excess may be present in the RGB LF of M13 (jCho et al.l 120051 ). These kinds of 
variations cannot be explained by the "global" physics that should apply to all globular 
cluster stars. These anomalies point toward fluctuations in the stellar initial mass function 
or composition-dependent effects. 



3.3.2. The RGB Bump 



A second feature of the LF presented here is a noticeable RGB bump. Typically, the 
RGB bump appears as a peak in the differential LF and as a change of slope in the cumulative 
LF (CLF). The bump provides a measure of the maximum depth reached by the outer 
convection zone during first dredge-up since it is the result of a pause in the star's evolution 
when the shell fusion source begins consuming material of constant, lower helium content 
( IFusi Pecci et al.lll990l ). Unfortunately, the number of stars occupying the bump gets smaller 
and the luminosity of the bump increases as the metallicity of the cluster decreases, making 
the bump harder to detect in metal-poor clusters. A small peak appears in our differential 
LF at V ~ 16.2, and a significant (2.5 — 3.5<r) change in slope occurs at the same position 
as the peak in the differential LF, as shown in Fig. [T71 

The relative brightness of the bump can be measured by comparing to th e ^-magnitude 



of the HB at the level of the RR Lyrae instability strip AVj^p p = V^^—Vhb ( Ferraro et al. 



19991 ). This indicator is a function of the total metallicity and the age of the cluster: an 
increase in m etallicity and /or a decrease in age are accompanied by a decrease in luminosity 
of the bump (IFerraro fc Montegriffoll2000l ). We find V& ump = 16.20 ± 0.05 mag, and Vhb = 
16.52 ± 0.11 (from interpolation between the average magnitudes of non-variable HB stars 
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at the blue and red e nds of the RR Lyrae i nstability strip), giving 

Av bum P = _o.32±o.l2. 

In the compilation of iFerraro et al. a zero-age HB reference point was calc ulated 

using the relation V ZAHB = V HB + 0.106[Fe/H] 2 + 0.236[Fe/H] + 0.193. IFerraro et al~1 found 
Vzahb = 16.62 ± 0.10, which is consistent with the value obtained here {Vzahb = 16.65 ± 
0.11). Our value of &Vzahb = — 0-45 ± 0.12 is considerably lower than tabulated values 
for other clusters with similar metallicities (M68: -0.60 ± 0.07 ; M92: -0.65 ± 0.12; M15: 



—0.65 ± 0.09). In a separate compilation, IZoccali et al.l (119991 ) measured a smaller value 
Av zahb = -°- 45 ± °- n for M 15, in better a greement with the value for NGC 5466. (The 
difference is primarily beca use IZoccali et al.l measured the bump position to be 0.16 mag 
fainter than IFerraro et al. I .) Clearly there is still some need for more precise comparisons of 
bumps in metal-poors clusters with theory. 

We believe, however, that the brightness of the bump should ultimately be judged using 
hydrogen-fusing stars as references because it avoids any effects of the poorly-understood 
physical processes (such as the helium flash and/or mass loss) associated with the creation 
of an HB star. Using the cluster LF (as seen in Figure IT41) . we again find that the observed 
RGB bump is fainter than model val ues by at least 0.3 mag when the models are shifted 
to match the cluster's main sequence. lHargis et al.l (120041 ) did similar comparisons between 
theoretical models and luminosity functions for M3, M5, M12, and M30 taken from the 
literature. With the exception of M30 (because the bump could not be identified), the 
position of the bump relative to the turnoff region agreed well with theory. NGC 5466 is 
thus the most metal-poor cluster this comparison has been done for. So at present we are left 
with the question of whether this might result from the cluster's low metallicity, or whether 
we have been the unfortunate victims of a fluctuation in the number of giant stars in this 
low-mass cluster. We therefore encourage the examination of the luminosity function of more 
massive metal-poor clusters to settle the question. 



3.3.3. Mass Function Exponent 



Two recent papers (IBelokurov et al.ll2006l ; iGrillmair fc Johnson! 120061 ) reported the dis- 
covery of ti dal streams covering many degrees around NGC 5466 in Sloan Digital Sky Sur- 



vey images. iGnedin fc Ostrikerl (119971 ) examined the Milky Way globular clusters and found 



NGC 5466 to be a cluster that has probably been strongly affected by disk shocking in the 
recent past. In our examination of the LF, we found a rather low value for the global main se- 
quence mass function slope. The mass function for a cluster is typically expressed as a power 
law (N(M) oc M~( 1+x ^), where the slope x = 1.35 is the standard Salpeter value. Generally, 
the present-day power-law index x varies from cluster to cluster. The mass function slopes 
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that best fit the upper LF of NGC 5466 around and above the MSTO have — 1 < x < 0. 
(Note that the best fit slope does depend on the models being used: the Yonsei-Yale models 
require a flatter slope than the Victoria-Regina and Teramo models.) Such a shallow mass 
function slope is unusual for a metal-poor cluster. For example, the cluster NGC 5053 has 
similar metallicity, position relative to the g alactic center and pla n e, and density str u cture , 
but still has a steep x ~ 2 mass function (IFahlman et al.lll99ll ). IDjorgovski et al.l (119931 ) 
found that mass function slopes in the range of 0.5 < < 0.8 are influenced primarily by 
the cluster's position in the galaxy, and to some extent by cluster metallicity. Based on both 
of those factors, NGC 5 466 should have a large r mass function slope (x ~ 3 according to the 
multivariate formula in IDjorgovski et al.lll993l ). Like other halo clusters, NGC 5466's orbit 
i s quite eccentric and will take the cluster more than 30 kpc away from the Galactic center 
(IDinescu et al.lll999l ). but it is currently on its way back into the halo after two relatively 
recent passes through the Galactic disk. Recent losses of low-mass stars may exp lain the 
recent identification of strong tidal tails near the cluster by lBelokurov et al.l (120061 ). 



4. Blue Stragglers 



Blue stragglers (BSSs) were first identified by ISandagd (119531 ) in the globular cluster 
M3. These stars are more massive than the turnoff mass and occupy the space in the CMD 
just bluer and brighter than the MSTO. Blue stragglers are found in clusters, and relatively 



more frequently in lower- luminosity clusters (IFerraro et al.l Il993l ; iPreston fc Snedenl 12000 



Piotto et al. 1120041 ; ISandquistl 120051 ) . From the various models proposed for the formation 
mechanism of BSSs, the "collision" theory (involving strong gravitational interactions be- 
tween previously unassociated single or binary stars) and the "mass-transfer" theory (in 
which the more massive star in a binary evolves and during its expansion transfers mass to 
its companion) are the strongest possibilities. There is a continuing interest in the study of 
BSSs because they may provide insight into the recent dynamical history of a cluster. 

In order to identify BSSs over the entire observed area of the cluster out to a radius 
of 11'.6, we used photometry from three datasets. In the core of the cluster we used the 
CFHT data p resented here for the first time. Outside of the CFHT field, we used the BV 
photometry of iJeon et al.ll2004l . which covered a field 1 1^6 on a side centered on the cluster. 
Finally, we used our KPNO data for the least crowded outskirts of the cluster. Even though 
NGC 5466 is a very low density cluster, the spatial resolution of the KPNO data was such 
that blends of stars would have resulted in the spurious identification of 10 objects as BSS 
in the intermediate portion of the field. 



48 BSS candidates were identified in NGC 5466 by 



Nemec &: Harris 



(119871 ). all located 
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between O'.l and 5f6 from the cluster center. In spite of the low cluster density, we find new 
BSS candidates at all radii and luminosity levels, and find several of thei r candidates are 
spurious. According to the CFHT photometry, the object with ID 45 from lNemec fc Harris 
is a blend of several fainter stars, no ne of which is a B SS. In addition, IDs 6 and 24 were 
identified as blends of stars using the IJeon et al.1 (12004 ) dataset. Our BSS list is presented 



m 
49^ 



able El The list includ es the nine known SX Phoenicis stars (ID 27, 29, 35, 38, 39, and 



Nemec & Mated Il990l : ID 3 (SX Phe 3), 36 (SX Phe 2), and 50 (SX Ph e l UJeon et al. 



20041 ) and the three eclipsing binaries (ID 19, 30, and 31 



Mateo et al. 



1990). New straggler 



candidates were given ID numbers that build upon the Nemec fc Harris! ( 1987 ) list. Fig 



mre 



[T8l shows the CMDs used to select the 94 identified BSSs in each of the three datasets. 



In order to use the BSSs to constrain cluster dynamics, we compared the normalized 
cumulative radial distribution of th e BSSs to the population of the giant branch, as shown 
in Fig. [19j iNemec fc Harris! (119871 ) found a 97.8% probability that their BSS sample was 



more centrally concentrated than red giants in the same magnitude range. Because their 
photometr y was taken in conditions of poorer seeing (compared to our CFHT photometry 
and that of lJeon et all ), their samples are likely to be somewhat incomplete near the cluster 
center. 

Our RGB sample contains 350 stars with magnitude V < 18.5. Kolmogorov-Smirnov 
(K— S) probability tests were used to test the hypothesis that both populations were drawn 
from the same parent population. The K— S probability that the BSSs are drawn from 
the same radial distribution as the RGB population is 8.1 x 10~ 7 , and 2.4 x 10 -4 for the 
comparison with the HB population. By contrast there is a probability of 0.27 that the RGB 
and HB samples are drawn from the same population. The concentration of BSSs toward 
the cluster center as compared to the RGB and HB samples is consistent with the idea that 
they are more massive than individual RGB stars, and as a result have been segregated by 
mass deeper within the cluster potential well. 



Piotto et al. I (120041 ) recently used samples of stragglers from the cores of 56 globular 
clusters to show that there was a strong correlation between F^g S = Nrgs/Nh b and inte- 



grated cluster V magnitude, and a weaker anti-correlation with central density. ISandquist 



(120051 ) examined an additional 13 low-luminosity globular clusters using similar selection 
criteria. NGC 5466 is an interesting cluster in relation to these samples because it has an 
integrated lu minosity that puts it at the faint end of the Piotto et al. sample (My t = —6.96; 
Harris! I1996T ). but with a central density that is nearly an order of magnitude lower than 



any of their clusters [log(/9 /(£v,eP c ~ 3 )) = 0.88] but comparable to clusters in the Sandquist 
sample. To put NGC 5466 in the context of the Piotto et al. and Sandquist samples, we 
selected a subset of our BSSs that satisfied the selection criteria in those studies (brighter 



-17- 



than the MSTO, and bl uer than th e MST O by 0.05 in B — V color). From mode fitting to 
the turnoff region in the I Jeon et all (120041 ) and CFHT data, we find (B — V)to — 0.367 and 
(V — I)r n = 0.511. A color offset of 0.05 in B — V corresponds to an offset of about 0.075 
in V - I JVandenBerg fc ClerrJ l2~0~0~3l ) for NGC 5466's metallicity. We find that 75 BSSs are 
brighter than the cluster turnoff (Vro — 19.99 ±0.05) and 0.05 bluer than the cluster turnoff 
in B — V. We have identified 97 HB stars in our datasets for NGC 5466, which gives a 
specific frequency F^ss = 0.77 ± 0.12 (with the error estimate from Poisson statistics). 



When compared to the iPiotto et al 



values (see Fig. ED), NGC 5466 falls within the 
general trend versus My t despite the cluster's low central density. On the other hand, NGC 
5466 has a lower Fgg S value than other cluste rs of similar central density (but lower total 
luminosity). As discussed by ISandquistl (120051 ). this provides additional evidence that the 
plateau in F^g S seen for clusters with log po ^ 2.5 is a result of the correlation between cluster 
integrated magnitude and central density. The lowest luminosity clusters in the Sandquist 
sample (E3 and Palomar 13) have central densities comparable to that of NGC 5466, but 
BSS frequencies that are several times higher. Another moderate-luminosity, low-density 
cluster (NGC 5053; Hiner et al., in preparation) similar to NGC 5466 shows a comparably 
low straggler frequency. BSSs produced via purely collisional means are no t likel y to show 
this kind of behavior. More likely is the scenario proposed by iDavies et al.1 (120041 ) in which 
binary stars that would normally produce BSSs are destroyed earlier in the cluster's history. 
More direct observational support for that hypothesis is needed though — for example, a 
detailed study of the variation of the binary star fractions as a function of integrated cluster 
magnitude. 

Fig. [21] shows that the number and frequency 



R 



BSS 



Nbss/Nbss 

^sample I ^sample 



of BSSs relative to the integrated V^-band flux (derived from a King model profile) as a func- 
tion of radius. Both frequencies increase toward the cluster center, and neith er shows signs 
of ris ing to ward larger radii. A s recent studies of denser c lusters show (M3: IFerraro et al. 
199.4 M55: IZaggia et al.lll997t 47 Tuc: IFerraro et allbooi NGC 6752: ISabbi et alil2004h . 
the BSS frequency generally decreases at i ntermediate radii and rises again at larger dis- 
tances. However, the cluster Palomar 13 (IClark et al.ll2004l ). which has a central density 
similar to NGC 5466, shows no sign of an increase in straggler frequency at large distance. 

In more massive clusters, the minimum in the BSS frequency is re ached approximatel y 
where the timescale for dynamical friction equals the age of the cluster (IWarren et al.ll2006l ). 
A similar calculation for the current structure of NGC 5466 indicates that this occurs at about 
270"(about 2.8r c ). This appears to be in the outer reaches of the core straggler distribution. 
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Because it is likely that NGC 5466 has lost a significant fraction of its mass, we expect 
that the current density structure of the cluster has not existed throughout its history and 
that NGC 5466 might have been able to dynamically relax stragglers to its core from larger 
distances earlier in its history. This may be showing that the global BSS population differs 
significantly between low-density/low-mass clusters and high-density/high-mass clusters. A 
lack of stragglers at large distance may be signature of large-scale tidal stripping of the 
cluster, which would remove both stragglers that would normally have formed in primordial 
binaries in the outer reaches of the cluster and ones that formed in the core but were given 
velocity kicks into orbits that would take them into the outer reaches. 

The case against a rise at large radius in Palomar 13 is stronger because the cluster 
has been surveyed out to 19 core radii, while in NGC 5466, we have only surveyed out to 
about 10 core radii (or 7.5 half-mass radii). Still, NGC 5466 probably should have an even 
more concentrated distribution of stragglers if its current density structure has existed for 
most of its history. In more massive clusters, the secondary rise in straggler frequency is 
observed between 8 and 10 r c . Unfortunately, further study of the stragglers in NGC 5466 
will probably be complicated by the strong tidal tails observed in the cluster. 

5. Conclusions 

Examinations of the luminosity functions of globular clusters continue to produce inter- 
esting tests of astrophysics. In this study, we found that NGC 5466 has a luminosity function 
that is in better overall agreement with theoretical models than the anomalous cluster M30, 
which has a similar low metallicity. In addition, we found that the relative numbers of red 
giant and main sequence stars may produce a fairly sensitive test of the physics near the 
core of red giants — specifically, helium diffusion and Coulomb interactions. However, we 
are not yet able to fully explain the differences between sets of theoretical models. 

Recent discoveries of large tidal tails associated with NGC 5466 suggest that this cluster 
has been strongly disrupted by interactions with the Galaxy. Our measured flat (—1 < x < 0) 
main-sequence luminosity function is unusual for a low-metallicity halo cluster. It is, however, 
consistent with the emerging picture of mass- segregation followed by tidal stripping. 

We have thoroughly re-examined the blue straggler population in the cluster, and de- 
tected a total of 94. The radial distribution of stragglers is clearly more centrally concen- 
trated than the RGB and HB populations. The frequency of blue stragglers in the cluster is 
relatively low — consistent with the observed anti-correlation between frequency and cluster 
luminosity, in spite of the cluster's very low central density. 



-19- 



We would like to thank the anonymous referee for helpful comments on the manuscript, 
Y. Jeon for providing us with an electronic copy of his photometric dataset, Y.-C. Kim for 
information on the Yonsei-Yale isochrones, and S. Cassisi for providing us with access to 
the Teramo set of models. This work has been funded through grants AST 00-98696 and 
05-07785 from the National Science Foundation to E.L.S. and M.B. 

REFERENCES 

Adelberger, E. G., et al. 1998, Rev. Mod. Phys., 70, 1265 

Bahcall, J. N., Pinsonneault, M. H., & Wasserburg, G. J. 1995, Rev. Mod. Phys., 67, 781 
Behr, B. B. 2003, ApJS, 149, 67 

Belokurov, V., Evans, N. W., Irwin, M. J., Hewett, P. C, & Wilkinson, M. I. 2006, ApJ, 
637, L29 

Bergbusch, P. A. 1996, AJ, 112, 1061 

Bjork, S. R., & Chaboyer, B. 2006, ApJ, 641, 1102 

Bolte, M. 1994, ApJ, 431, 223 

Bono, G., Cassisi, S., Zoccali, M., & Piotto, G. 2001, ApJ, 546, L109 
Buries, S., Nollett, K. M., & Turner, M. S. 2001, ApJ, 552, LI 
Carretta, E., & Gratton, R. G. 1997, A&AS, 121, 95 

Cassisi, S., Salaris, M., Castelli, F., & Pietrinferni, A. 2004, ApJ, 616, 498 
Cassisi, S., Salaris, M., & Irwin, A. W. 2003, ApJ, 588, 862 
Cho, D.-H., Lee, S.-G., Jeon, Y.-B., & Sim, K. J. 2005, AJ, 129, 1922 
Clark, L. L., Sandquist, E. L., & Bolte, M. 2004, AJ, 128, 3019 

D'Antona, F., Bellazzini, M., Caloi, V., Pecci, F. F., Galleti, S., & Rood, R. T. 2005, ApJ, 
631, 868 

Davies, M. B., Piotto, G., k de Angeli, F. 2004, MNRAS, 349, 129 
degl'Innocenti, S., Weiss, A., & Leone, L. 1997, A&A, 319, 487 



-20- 

Demarque, P., Woo, J.-H., Kim, Y.-C, & Yi, S. K. 2004, ApJS, 155, 667 
Dinescu, D. L, Girard, T. M., & van Altena, W. F. 1999, AJ, 117, 1792 
Djorgovski, S., Piotto, G., & Capaccioli, M. 1993, AJ, 105, 2148 
Eggleton, P. P., Faulkner, J., & Flannery, B. P. 1973, A&A, 23, 325 
Fahlman, G. G., Richer, H. B., & Nemec, J. 1991, ApJ, 380, 124 

Ferraro, F. R., Beccari, G., Rood, R. T., Bellazzini, M., Sills, A., & Sabbi, E. 2004, ApJ, 
603, 127 

Ferraro, F. R., Messineo, M., Fusi Pecci, F., de Palo, M. A., Straniero, O., Chieffi, A., & 
Limongi, M. 1999, AJ, 118, 1738 

Ferraro, F. R. & Montegriffo, P. 2000, A J, 119, 1282 

Ferraro, F. R., Paltrinieri, B., & Cacciari, C 1999, Mem. Soc. Astron. Italiana, 70, 599 

Ferraro, F. R., Fusi Pecci, F., Cacciari, C, Corsi, C, Buonanno, R., Fahlman, G. G., & 
Richer, H. B. 1993, AJ, 106, 2324 

Fusi Pecci, F., Ferraro, F. R., Crocker, D. A., Rood, R. T., & Buonanno, R. 1990, A&A, 
238, 95 

Gallart, C, Zoccali, M., & Aparicio, A. 2005, ARA&A, 43, 387 
Gnedin, O. Y., & Ostriker, J. P. 1997, ApJ, 474, 223 
Grillmair, C. J., & Johnson, R. 2006, ApJ, 639, L17 

Guenther, D. B., Demarque, P., Kim, Y.-C, & Pinsonneault, M. H. 1992, ApJ, 387, 372 

Guhathakurta, P., Webster, Z. T., Yanny, B., Schneider, D. P., & Bahcall, J. N. 1998, AJ, 
116, 1757 

Hargis, J. R., Sandquist, E. L., & Bolte, M. 2004, ApJ, 608, 243 

Harpaz, A., & Kovetz, A. 1988, ApJ, 331, 898 

Harris, W. E. 1996, AJ, 112, 1487 

Hathaway, D. H. 1996, ApJ, 460, 1027 

Jeon, Y., Lee, M. G., Kim, S., & Lee, H. 2004, AJ, 128, 287 



- 21 - 



Johnson, J. A., & Bolte, M. 1998, AJ, 115, 693 

Kim, Y., Demarque, P., Yi, S. K., & Alexander, D. R. 2002, ApJS, 143, 499 
Landolt, A. U. 1992, AJ, 104, 340 

Lejeune, T., Cuisinier, F., & Buser, R. 1998, A&AS, 130, 65 

Mateo, M., Harris, H. C, Nemec, J., & Olszewski, E. W. 1990, AJ, 100, 469 

McCarthy, J. K., & Nemec, J. M. 1997, ApJ, 482, 203 

Nemec, J. M. & Harris, H. C. 1987, ApJ, 316, 172 

Nemec, J., & Mateo, M. 1990, ASP Conf. Ser. 11: Confrontation Between Stellar Pulsation 
and Evolution, 11, 64 

Norris, J. E. 2004, ApJ, 612, L25 

Olive, K. A., & Skillman, E. D. 2004, ApJ, 617, 29 

Piotto, C, De Angeli, F., King I. R., Djorgovski, S. C, Bono, C, Cassisi, S. , Meylan, C, 
Recio-Blanco, A. , Rich, R. M. & Davies, M. B. 2004, ApJ, 604, 109 

Piotto, C, & Zoccali, M. 1999, A&A, 345, 485 

Pollard, D. L., Sandquist, E. L., Hargis, J. R., & Bolte, M. 2005, ApJ, 628, 729 

Preston, G. W., & Sneden, C. 2000, AJ, 120, 1014 

Pritzl, B. J., Venn, K. A., & Irwin, M. 2005, A J, 130, 2140 

Proffitt, C. R., & Michaud, G. 1991, ApJ, 371, 584 

Proffitt, C. R., & Vandenberg, D. A. 1991, ApJS, 77, 473 

Pryor C. , McClure, R. D., Fletcher, J. M., & Hesser, J. E., 1991, AJ, 102, 1026 

Ratcliff, S. J. 1987, ApJ, 318, 196 

Renzini, A. & Fusi Pecci, F. 1988, ARA&A, 26, 199 

Rogers, F. J., & Nayfonov, A. 2002, ApJ, 576, 1064 

Rogers, F. J., Swenson, F. J., & Iglesias, C. A. 1996, ApJ, 456, 902 

Rood, R. T., et al. 1999, ApJ, 523, 752 



-22 - 

Rosenberg, A., Aparicio, A., Saviane, L, & Piotto, G. 2000, A&AS, 145, 451 
Sabbi, E., Ferraro, F. R., Sills, A., & Rood, R. T. 2004, ApJ, 617, 1296 
Salaris, M., Riello, M., Cassisi, S., & Piotto, G. 2004, A&A, 420, 911 
Salaris, M., & Weiss, A. 2002, A&A, 388, 492 
Sandage, A. R. 1953, AJ, 58, 61 

Sandquist, E. L., Bolte, M., Langer, G. E., Hesser, J. E., & Mendes de Oliveira, C. 1999, 
ApJ, 158, 262 

Sandquist, E. L., Bolte, M., Stetson, P. B., & Hesser, J. E. 1996, ApJ, 470, 910 
Sandquist, E. L. 2005, ApJ, 635, 73 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 
Searle, L., & Zinn, R. 1978, ApJ, 225, 357 
Spergel, D. N., et al. 2006, submitted 
Stetson, P. B. 2000, PASP, 112, 925 

Stetson, P. B. 1991, ASP Conf. Ser. 13: The Formation and Evolution of Star Clusters, 13, 
88 

Stetson, P. B. 1990, PASP, 102, 932 
Stetson, P. B. 1987, PASP, 99, 191 

Straniero, O., Chieffi, A., & Limongi, M. 1997, ApJ, 490, 425 

VandenBerg, D. A., Bergbusch, P. A., & Dowler, P. D. 2006, ApJS, 162, 375 

Vandenberg, D. A., Bolte, M., & Stetson, P. B. 1990, AJ, 100, 445 

VandenBerg, D. A., & Clem, J. L. 2003, AJ, 126, 778 

Vandenberg, D. A., Larson, A. M., & de Propris, R. 1998, PASP, 110, 98 

VandenBerg, D. A., Richard, O., Michaud, G., & Richer, J. 2002, ApJ, 571, 487 

VandenBerg, D. A., Swenson, F. J., Rogers, F. J., Iglesias, C. A., & Alexander, D. R. 2000, 
ApJ, 532, 430 



-23- 

Warren, S. R., Sandquist, E. L., & Bolte, M. 2006, ApJ, 648, 1026 

Zaggia, S. R., Piotto, G., & Capaccioli, M. 1997, A&A, 327, 1004 

Zoccali, M., Cassisi, S., Piotto, G., Bono, G., & Salaris, M. 1999, ApJ, 518, L49 

Zoccali, M., & Piotto, G. 2000, A&A, 358, 943 

Zinn, R. & West, M. J. 1984, ApJS, 55, 45 

Zinn, R. 1980, ApJ, 241, 602 



This preprint was prepared with the AAS IATgX macros v5.2. 



-24- 



0.1 
0.05 

< u 
-0.05 

-0.1 



0.1 
0.05 




L 0.006 (0.027) 

= * xx *x x **» ** U 


, 1 1 1 , x! ^ 

XX X X 

< X X X X 

< s, x x x JXx 


i i i i i i i * i i i i i x i i i i i 

X % x x 

x S V x 

-X X M *&<X*3* X x* 


IX xx ^ * *x x x 

-1 1 1 1 1 1 


X x 

x x xx x 

XX — 

I ill r 


x>** x #ir«x 

x-x x xx x 


12 14 

B 


16 18 0.5 1 1.5 

B-V 



-0.05 -- 
-0.1 



0.1 
0.05 
< 
-0.05 f 
-0.1 



7 o.i 

> 



- 1 1 1 1 1 1 1 1 1 1 1 - 

L 0.001 (0.016) xj 

=-« * * x ^xyx x xx | x ^ x x x gxx >; : 


i i i i i i i i i i i i i i i i 

~r x * « » x 


; x x x fx x ^ jj x x x ^ ^ ^ ^ _ 
- 1 1 1 1 1 1 1 1 1 1 1 1 1 - 


_ ^ X X 


12 14 16 18 0.5 1 1.5 

V B-V 



<3 
CQ 



0.05 


0.05 

-0.1 



- 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -1 

L 0.005 (0.026) x x x x x x x X J 

: „ x x x x fe« x« X ! * x : 

»» X X X X g X * * X X x X * 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

* X 
X X X 

x y£ x x x - 

XX^X X X 


*x |!S xj *x x 

-i " 8 x x : 

~ X ~ 

- i i i 1 i i r 1 i i i 1 i i - 


-^r^ &P*feP x a g — * 

* x J x 

_ X _ 

[ 1 I I 1 1 1 1 [ 1 1 1 I 1 1 1 1 1 1 1 1 1 1 I 1 


12 14 16 0.5 1 1.5 2 

I V-I 


_ 1 1 1 | 1 1 1 1 1 1 1 1 | 1 1 1 1 1 1 1 1 _ 

L 0.007 (0.02a) x J 

= X xxXj| xx x X * X z 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

: -0.004 (0.021) x 7 

x x ** 


X **« X * X X 
" 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 I 1 1 1 1 " 


: «|; ! xfxHx x x x x *; 

X x *x x xX 



0.5 1 

(B-V) 



1.5 



0.5 1 

(v-i) 



1.5 



Fig. 1.— Pho tometric residuals (in the sense of this study minus those of lLandoftlll992l and 



Stetson! 120001 ) of primary standard stars. The median residuals are listed in the panels with 



the semi-interquartile range (half the magnitude difference between the 25% and 75% points 
in the ordered list of residuals) given in parentheses. 
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Fig. 2. — Photometric residuals (in the sense of the final PSF photometry minus standard 
aperture photometry values) of secondary standard stars. 
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Fig. 3. — Residuals (in the sense of this study minus IStetsonl l2000l ) from the star-by-star 
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Fig. 4. — Residuals (in the sense of this study minus iRosenberg et al.ll2000l ) from the star- 
by-star comparison. The median residuals and the plots versus color have been restricted 
to brighter stars (V < 20 and / < 19) to make the comparisons clearer. The numbers in 
parentheses are the semi-interquartile ranges (see FigJT]). 
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Fig. 5. — Residuals (in the sense of this study minus iJeon et al.ll2004l ) from the star-by-star 
comparison. The median residuals and the plots versus color have been restricted to brighter 
stars (B < 20 and V < 19.5) to make the comparisons clearer. The numbers in parentheses 
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Fig. 6. — External V magnitude errors cr ext (V) as a function of radius and magnitude 
determined from artificial star tests, with exponential fits shown by the solid lines. 
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Fig. 7. — Magnitude biases S(V) determined from artificial star tests as a function of radius 
and magnitude. 
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Fig. 8. — Total recovery probability F(V) determined from artificial star tests as a function 
of radius and magnitude. 
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Fig. 9. — Completeness fraction f(V) determined from artificial star tests as a function of 
radius and magnitude. 
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Fig. 10. — Color-magnitude diagrams for all stars measured in the KPNO and CFHT images. 
The BV fiducial (Table [3]) is also plotted in the left panel. 
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Fig. 11. — Comparison of the observed fiducial sequence of NGC 5466 with the isochrones 
of the Teramo, Victoria-Regina, and Yonsei- Yale groups. The isochrones have been shifted 
horizontally so that the turnoff colors align, and shifted vertically to align the main sequence 
point 0.05 mag redder than the turnoff. On the giant branch, the ages increase from the 
reddest to the bluest isochrone. 
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Fig. 12. — Same as Fig. [TT], except for more metal- rich models. 
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Fig. 13. — Comparison of the observed V— band LF of NGC 5466 with theoretical models 
of the Yonsei-Yale, Teramo, and Victoria- Regina groups assuming (m — M) v = 16. 
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Fig. 14. — Comparison of the observed V— band LF of NGC 5466 with theoretical models 
of the Victoria-Regina, Yonsei-Yale, and Teramo groups using magnitude shifts that bring 
the main sequence point 0.05 mag redder than the turnoff into alignment. The models have 
been normalized to the two bins on either side of the turnoff (V = 19.94). 
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Fig. 15. — Comparison of the observed LFs of M3, M5, M12, and M30 with theoretical 
models of the Teramo group using magnitude shifts that bring the main sequence point 0.05 
mag redder than the turnoff into alignment. The models have been normalized to bins on 
either side of the turnoff. 
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Fig. 16. — Fractional difference between the observed RGB-MS number ratios (for M5, M12 
M3, M30, and NGC 5466, from left to right) and the theoretical predicti ons from the Yonsei- 
Yale, Teramo, and Victoria- Regina models. The le ft panels use the 



Carretta fc Gratton 



(119971 ) metallicity scale, and the right panels use the IZinn fc Westl ( 119841 ) scale. The sense 
is (observed — theoretical) / observed. 
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Fig. 17. — The cumulative luminosity function for bright RGB stars derived from the pho- 
tometry of Jeon et al. (2004; 286 stars) and from the KPNO dataset (338 stars) presented 
here. Dotted lines show fits to the data for stars above and below the position of the apparent 
bump. 
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Fig. 18. — Blue straggler selection for NGC 546 6. The stars plotte d in each panel show the 
entire sample used for the selection: stars from iJeon et al.l (120041 ) in the middle panel are 
only those s tars outside the CFHT field, and KPNO stars in the right pan el are only those 
outsid e the IJeon et al.l field. Open squares show stragglers identified by iNemec fc Harris 
(119871 ). and open circles are new candidates. 
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Fig. 19. — Normalized cumulative radial distributions for RGB stars (dashed line), HB stars 
(dotted line), and BSSs (solid line). 
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Fig. 20. — Relative frequencies of blue stragglers as a function of cluster absolute magnitude 
and central den s ity. T he solid square is NGC 5466, the open squares ar e globular clusters 
from Sandauist (2005), and all other noints are from Piotto et al. (2004). Open circles are 
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Fig. 21. — Frequency of BSS relative to the integrated V^-band flux of detected cluster stars 
(top panel) and specific frequency of blue stragglers relative as a function of radius (bottom 
panel). 
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Table 1. Photometric Observation Log for NGC 5466 



UT Date 


Filters 


N 


Exposure Time (s) 


Airmass 


1995 May 4 


By 


1,1 


60 


1.01,1.12 


1995 May 4 


B 


2 


300 


1.03, 1.03 


1995 May 4 


By 


2,1 


600 


1.0,1.11,1.0 


1995 May 5 


B,I 


2,2 


300 


1.03,1.01,1.02,1.01 


1995 May 9 


B,V,I 


1,1,1 


60 


1.12,1.13,1.14 



Table 2. V-Band Luminosity Function 



V logiV G hi gh &low 



13 


.532 





.5005 





.1761 


0. 


.3010 


13 


.832 


-0. 


.1015 





.3010 
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.0000 


14. 


.132 


-0. 


.1015 





.3010 
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.0000 


14. 


.432 





.5975 





.1605 





.2575 


14. 


.732 





.8985 





.1193 


0. 


.1651 


15. 


,032 





.6767 


0. 


.1487 





,2279 


15. 


,332 





,6767 





.1487 
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.2279 


15. 


.632 





.5976 





.1606 





,2575 


15. 


,932 
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,1651 


16 


.232 
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.2410 
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.1041 


16 


.532 


1 
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.832 
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.6996 
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.0201 


20 


.132 
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.7738 





.0178 
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.432 
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21 


,032 
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.0036 
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,0146 


21 


,332 
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,0496 





.0137 
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21 
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0. 


.0132 





.0136 
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Table 3. Fiducial sequence for NGC 5466 
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0. 


636 
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17 


.820 


0. 


.647 


31 


17. 


.586 


0. 


.654 


23 


17. 


.403 


0. 


.672 


20 


17. 


.194 


0. 


.677 


19 


16 


.988 


0. 


693 


11 


16 


.784 


0. 


.722 


10 


16 


.613 


0. 


.725 


6 


16 


.400 


0. 


.752 


12 


16 


.199 


0. 


.771 


16 
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Table 3 — Continued 



V 


B-V 


iV a 


16.066 


0.779 


5 


15.825 


0.809 


7 


15.602 


0.827 


4 


15.407 


0.857 


6 


15.111 


0.917 


1 


15.006 


0.930 


5 


14.786 


0.946 


2 


14.590 


0.986 


11 


14.438 


1.044 


2 



a Number of stars 
used to determine 
fiducial point 
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Table 4. RGB-MSTO Number Ratios 



Source a TO Sample RGB Sample Y N rgb /N M sto [Fe/H] 



NGC 5466 


19.682 < V < 20.282 16.982 < V < 18.482 




0.162 ±0.013 




VR 




0.235 


0.132 


-2.22 


VR 




0.235 


0.130 


-2.14 


T 




0.245 


0.148 


-2.22 


T 




0.245 


0.146 


-2.14 


YY 




0.230 


0.162 


-2.22 


YY 




0.230 


0.160 


-2.14 


M3 


18.80 < V < 19.40 16.40 < V < 18.00 




0.168 ±0.008 




VR 




0.235 


0.170 


-1.66 


VR 




0.235 


0.158 


-1.34 


T 




0.246 


0.181 


-1.66 


T 




0.248 


0.170 


-1.34 


YY 




0.230 


0.185 


-1.66 


YY 




0.230 


0.162 


-1.34 


M5 


19.13 < B < 19.73 16.33 < B < 17.93 




0.110 ±0.006 




VR 




0.235 


0.106 


-1.40 


VR 




0.235 


0.097 


-1.11 


T 




0.248 


0.114 


-1.40 


T 




0.251 


0.106 


-1.11 


YY 




0.230 


0.119 


-1.40 


YY 




0.230 


0.105 


-1.11 


M12 


18.14 < V < 18.74 15.59 < V < 17.24 




0.158 ±0.011 




VR 




0.235 


0.144 


-1.40 


VR 




0.235 


0.118 


-1.14 


T 




0.248 


0.155 


-1.40 


T 




0.251 


0.144 


-1.14 


YY 




0.230 


0.150 


-1.40 


YY 




0.230 


0.137 


-1.14 


M30 


18.33 < V < 18.93 15.78 < V < 17 A3 




0.214 ±0.017 




VR 




0.235 


0.174 


-2.13 


VR 




0.235 


0.153 


-1.91 
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Table 4 — Continued 



Source a TO Sample 


RGB Sample Y 


Nrgb/Nmsto 


[Fe/H] 


T 


0.245 


0.173 


-2.13 


T 


0.246 


0.165 


-1.91 


YY 


0.230 


0.194 


-2.13 


YY 


0.230 


0.180 


-1.91 



a VR: Victoria- Regina models (no diffusion); YY: Yonsei-Yale models (He 
diffusion); T: Teramo models (no diffusion). All models are for an age of 
12 Gyr. 
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Table 5. Selected Star Populations in NGC 5466 



ID Aq(") A6(") B a B V o v I oi Alternate ID Ref. f 



Blue Stragglers 



1 


137 


74 


32 


86 


19 


1477 






19 


0038 














605 


J 


1 


137 


74 


32 


86 


19 


1446 





0122 


19 


0182 





0280 


18 


7606 





0400 


809 


K 


2 


-12 


18 


-5 


60 










18 


7591 





0116 


18 


5972 





0061 


10313 


C 


2 


-12 


18 


-5 


60 


18 


8461 





0147 


18 


7209 





0258 


18 


6245 





0504 


2176 


K 


3 


-8 


68 


14 


99 










19 


1710 





0100 


18 


9595 





0059 


9339 


C 


3 


-8 


68 


14 


99 


19 


4373 





0155 


19 


2441 





0296 


19 


1073 





0563 


2129 


K 


4 


-77 


31 


83 


56 


19 


2057 






19 


0546 














646 


J 


4 


-77 


31 


83 


56 


19 


2266 





0123 


19 


0548 





0256 


18 


8465 





0460 


2880 


K 


5 


-81 


70 


65 


06 


18 


5511 






18 


2697 














389 


J 


5 


-81 


70 


65 


06 


18 


5632 





0121 


18 


2740 





0243 


17 


8309 





0255 


2895 


K 


7 


-132 


96 


-3 


18 


18 


8010 






18 


7046 














504 


J 


7 


-132 


96 


-3 


18 


18 


8482 





0113 


18 


7379 





0217 


18 


5538 





0386 


3289 


K 


8 


-144 


85 


3 


07 


18 


8436 






18 


6640 














498 


J 


8 


-144 


85 


3 


07 


18 


8726 





0120 


18 


6345 





0234 


18 


3367 





0298 


3358 


K 


9 


-90 


90 


83 


38 


19 


1757 






19 


0256 














626 


J 


9 


-90 


90 


83 


38 


19 


1385 





0144 


18 


9711 





0257 


18 


7508 





0387 


2973 


K 


10 


-44 


94 


96 


04 


19 


3627 






19 


1882 














751 


J 


10 


-44 


94 


96 


04 


19 


4145 





0139 


19 


2190 





0261 


18 


9077 





0494 


2520 


K 



a Photometry Sources: K: KPNO data from this paper, C: CFHT data from this paper, J: Jeon et aD ([' 



Note. — The complete version of this table is in the electronic edition of the Journal. The printed edit: 
only a sample. 



